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We calculate the pion-pion elastic scattering phase shifts for pion masses from the chiral limit 
to values of interest for lattice studies. At low energies, we use the standard Chiral Perturbation 
Theory expressions to one and two loops. In addition, we study the phase shifts mass dependence 
in the resonance region by means of dispersion theory in the form of unitarized Chiral Perturbation 
Theory and the Inverse Amplitude Method. We pay particular attention to the case when resonances 
are close to threshold, illustrating the different behavior between scalar and vector resonances. We 
also provide the estimation of uncertainties, which are dominated by those of the O(p^) chiral 
parameters. 
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I. INTRODUCTION. 

Elastic pion-pion scattering has been an object of study 
for many decades due to several reasons. In particular, 
pions are very relevant in the description of final states 
in other hadronic processes. Also, the two pion corre- 
lated exchange in the scalar-isoscalar channel is the main 
contribution to nucleon-nucleon attraction, and has been 
interpreted for long as a scalar "sigma" resonance [H, 
whose existence, mass and width have been the subject 
of an intense debate. Actually this resonance, nowadays 
called /o(600), appears as a pole deep in the second Rie- 
mann sheet of the scattering amplitude (see the "Note 
on scalar mesons" in for a detailed account). Finally, 
the pion-pion interaction at low energies is also relevant 
for the determination of light quark mass ratios and the 
size of the chiral condensate 

On the theory side, unfortunately, neither the elastic 
resonance region nor the low energy region are accessible 
to perturbative QCD calculations. In order to describe 
these processes in terms of quarks and gluons one should 
rely on lattice techniques. For a long time, these tech- 
niques have found little applications in this low-energy 
realm due to complications on the implementation of 
chiral symmetry, the small physical values of the light 
quarks and other technicalities as the existence of quark- 
line disconnected diagrams in some channels. However, 
very recently, lattice results have become available for 
the /o(770) and /n(600) resonance masses the pion 

decay constant [3, [13, [lH or even the isospin 2 scattering 
length (ill . [l^ , obtained with pion masses which are not 
too far from the physical values. Recent developments 
[T^ l in algorithms may make disconnected diagrams for 
multi-hadron calculations tractable in the not too distant 
future. This means that pion-pion scattering phase shifts 
might be calculable soon within lattice QCD. Actually, 
some first results for the isospin 2 waves have been ob- 
tained for still somewhat large pion masses [13, [IE] • Of 
course, lattice calculations still have systematic uncer- 
tainties which are hard to estimate and they always rely 
on modified actions, finite volumes, and other complica- 



tions so that their physical results are actually extrapo- 
lations to the physical limit. It is therefore necessary to 
understand how these chiral or physical extrapolations 
should be carried out. 

Fortunately, even though we cannot rely on perturba- 
tive QCD at low energies, we can still use its effective 
low ener gy t heory, known as Chiral Perturbation Theory 
(ChPT) [l6|, which provides a rigorous, systematic and 
model independent expansion of hadronic observables in 
terms of the external meson momenta and the relatively 
small pion mass. We will very briefly review ChPT in 
section [m mostly to introduce the required notation. 

Within ChPT, the quark mass dependence appears in 
a model independent way through the pion mass squared, 
which is also described as an expansion. Remarkably, the 
isospin 1 — 2 scattering length dependence found on 
the lattice is rather well described by just leading order 
ChPT up to surprisingly large pion masses ^1, 12] and 
the one-loop corrections seem to be rather small. In this 
work we will first study the evolution of the lowest five 
pion-pion scattering phase shifts, with definite isospin 
and angular momentum (/, J) = (0, 0), (1, 1), (2, 0), (0, 2) 
and (2,2), using the one and two-loop standard ChPT 
expressions, estimating the uncertainties due to the rel- 
atively poor knowledge of the low energy constants — 
particularly those at two loops. Of course, this approach 
is limited to low masses and momenta and cannot be used 
to describe resonances, although, in principle it should 
be able to describe their low energy tails, through, for 
instance, the low energy scattering phase shifts. This is 
the reason why one of the aims of this work is to study the 
evolution of all tttt scattering phase shifts at low energy 
within standard ChPT. 

Beyond the low energy regime, it is still possible to ob- 
tain the quark mass dependence of hadronic observables, 
by combining ChPT with dispersion relations. Thus, 
in section |V| we briefly review the Inverse Amplitude 
Method (lAM) (TtI - Ii^ , obtained by using the elastic ap- 
proximation together with ChPT, to calculate the sub- 
traction constants and the left cut contribution of a dis- 
persion relation for the inverse of the partial waves. This 
technique provides a description of meson-meson scatter- 
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ing which is simultaneously compatible with the ChPT 
low energy description but also generates the lightest 
elastic resonance on each channel. By applying this tech- 
nique to the TTTT scattering amplitude to one-loop in SU(2) 
ChPT, some of us have calculated the pion mass depen- 
dence of the p(770) and /o(600) masses and widths [2l] |. 
Interestingly, this method had already been applied to 
study only the /o(600) quark mass dependence and its 
influence, through the nucleon-nucleon interaction, on 
the production of carbon and oxygen and its anthropic 
implications [l^l- Recently [1^, some of us have also 
calculated the k(800) and if* (892) mass and width de- 
pendence with respect to the non strange-quark mass, as 
well as the dependence of all these four resonances with 
respect to the strange quark mass. And even more re- 
cently [13] we have extended to too loops the analysis 
of the p(770) and /o(600) resonances within unitarized 
elastic tttt scattering. 

The lAM results for the dependence of the p(770) 
agree nicely with the estimations for the two first coeffi- 
cients of its chiral expansion [2^ , and also with the ex- 
isting lattice results [J-[8|. The comparison with lattice 
is relatively straightforward in this case since the yo(770) 
is not extremely wide and it is actually calculated as a 
state of the spectrum. 

Unfortunately, the comparison of the lAM with lat- 
tice results will not be so straightforward for the scalar 
channels. First, we find of particular interest the repul- 
sive 1 = 2 channels. Note that these channels have no 
resonances, so that neither the spectroscopic studies on 
the lattice nor our pole studies with the lAM [2ll.[2l[23| 
address this case. However, this is the simplest channel 
for scattering lattice studies and, as commented above, 
there are already some lattice results for the scattering 
length down to relatively low pion masses [ill [l2l] and 
for phase-shifts but only for TOtt — 400 MeV or higher 

Mm 

Second, we are also interested in the much debated 
isoscalar channel. Of course, given the status of the a 
or /o(600), reliable lattice results would be most wel- 
come. Unfortunately lattice calculations in this channel 
are hard due to disconnected diagrams, but also their 
interpretation would be complicated because this reso- 
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nance is extremely wide (see 
In addition, it was shown in (21|, 
high masses, the /o(600), being a scalar, becomes a vir- 
tual state - a pole in the second Riemann sheet below 
threshold - which is not a physical state of the spec- 
trum. Therefore, since spectroscopic (or "pole") lattice 
studies of the a may be rather complicated, a study of 
the scalar phase shift, as the one presented here, deserves 
more interest. 

These are the motivations to study the chiral extrapo- 
lation of phase shifts either from standard or unitarized 
ChPT. This will be done first for standard ChPT to next 
to leading order (NLO) in section UlI Al and then to next 
to next to leading order (NNLO) in section IIIIBI Sur- 
prisingly, in both cases, the predicted behavior for the 



phase shift in the p(770) may look counterintuitive when 
compared with present lattice calculations of the p(770) 
mass TO,r dependence. This discussion deserves a sepa- 
rated section, in which we also evaluate the pion mass 
dependence of the "size" of the p(770). Next we will 
present the lAM resuhs for NLO ChPT in section IVI Al 
and for NNLO in section lyTBl We will discuss and sum- 
marize all our findings in section rVIIII 



II. CHIRAL PERTURBATION THEORY 

Pions are the Goldstone bosons associated to the spon- 
taneous chiral symmetry breaking of QCD. If quarks were 
strictly massless, pions would be massless too and sep- 
arated by a gap of the order of 1 GeV from the rest of 
hadrons, becoming the relevant QCD low energy degrees 
of freedom. Chiral Perturbation Theory (ChPT) [Ig] is 
nothing but the most general Lagrangian built out as an 
expansion in pion momenta (i.e., derivatives) respecting 
the QCD symmetries. In real life, though, the u and d 
quarks have a very small mass, that we will take in the 
isospin limit as to = (to„ -I- TOd)/2, which can be treated 
as a perturbation within ChPT. As a consequence pi- 
ons have a physical mass of to^ = 139.57 MeV, whose 
model independent perturbative expansion in terms of 
TO is given by ChPT. In summary, the QCD low energy 
theory we will use is SU(2) ChPT [161, which corresponds 
to considering the u and d quarks only and integrating 
out the other four quarks, whose effect will be included 
in the low energy constants (LECs) that multiply each 
term of the ChPT Lagrangian. In this way only pions 
will circulate in the loops. Hence, by varying the pion 
mass while keeping the ChPT low energy constants fixed, 
we are sensitive to the light quark mass dependence for 
constant s, c, b and t masses. 



A. Perturbative nn scattering within ChPT 

Pion-pion elastic scattering is customarily described in 
terms of partial wave amplitudes i j^"* (s) of definite isospin 
/ and angular momentum J, where s is the Mandelstam 
variable, although for simplicity we will drop these in- 
dices when there is no possible confusion. From ChPT 
these partial waves are obtained as a series expansion 
t = t2 +t4 + tQ • • • , with tk = 0(p/47r/^)''', where p stands 
generically for center of mass momenta or pion masses. 
The leading order (LO) t2 is 0{p^) and is universal [265 '^^ 
the sense that it only depends on the scale /^r ~ 92.4 MeV 
and TOtt. The NLO calculation yields t4 ,163 ob- 
tained from one-loop diagrams with LO vertices and tree 
diagrams from the NLO Lagrangian terms, which are 
multiplied by some low energy constants (LECs), called 
Z[(^). These LECs absorb the dependence on the loop 
regularization scale fj,, and are determined by the un- 
derlying QCD dynamics. Their measured values can be 
found in Table HI Something similar happens with the 
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0(p*) LECs (xlO^^) 


0(p ) 


LJiOs(xlU j 




-3.98 ± 0.62 




-0.60 ± 0.35 




1.89 ± 0.23 




1.28 ± 0.74 




0.18 ± 1.11 




-1.68 ± 0.97 




6.17 ± 1.39 




-1.00 ± 0.58 






rlM 


1.52 ± 0.42 






reM 


0.40 ± 0.04 








0.00 ± 1.20 



TABLE I: ChPT low energy constants from [2g] that con- 
tribute to TVTT scattering to O(p^) and 0{p^) that we use in 
our standard ChPT calculations. The value for Isdj,) comes 
from a recent analysis of the lattice results [s^]. The renor- 
malization scale is set to fj, — 770 MeV. Errors are only statis- 
tical or "only account for the noise seen in the calculations" 
of [2^. The first four and their uncertainties are obtained 
from resonance saturation. The r^ifJ-) value is from [3^ . 



NNLO result te [27[, which has two loop contributions 
with LO vertices, one-loop contributions with one LO 
vertex and one NLO vertex containing some li, plus tree 
level diagrams with NNLO vertices, whose LECs appear 
only in six combinations now called r^ifJ-), whose esti- 
mated values are listed also in Table HI All these LECs 
carry a scale dependence that cancels that from loop inte- 
grals, so that observables are scale independent and finite 
order by order. 

Let us remark that we write the tttt scattering ampli- 
tude in terms of the physical constants and , which 
are obtained as expansions in powers of the LO pion 
mass. Actually, I3 and I4 appear at NLO in tttt scatter- 
ing through these ttXtt and expansions, but in contribu- 
tions that depend stronger on the pion mass and softer on 
the energy than those containng the other LECs. Thus I3 
and I4 are harder to determine experimentally and have 
the largest uncertainty. This is particularly severe for 
and that is why we have used its lattice determination 
^ quoted in Table H 

At NNLO, the expansion of f-^ on the physical pion 
mass requires an additional parameter r f , also listed 
in Table H] Note that there is an additional 0{p^) con- 
stant, rM, which appears in the NNLO chiral expansion 
of the physical pion mass TOtt in powers of the quark mass 
m, but such a constant would only be needed in order to 
study the quark mass dependence of observables. How- 
ever, quark masses carry some renormalization scale and 
scheme dependence and most lattice results provide their 
results in terms of the physical pion mass. That is why 
here we will study the dependence of scattering phases 
on the physical pion mass and not on the quark mass. 
Therefore we do not need tm- 

We show in Table |T] the estimated statistical uncertain- 
ties of the LECs (for , rg they are described as the noise 
in the dispersive calculation of (1^). Systematic uncer- 
tainties are large and harder to estimate; for illustration 
we also provide in Table |lT] other values found in the lit- 
erature at OIp^) and 0{p^). We consider the spread on 
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1 ci^ 1^ 
lU li 


1 n3 1^ 

LU I2 


iU <3 


iU t4 


/^l 1~>rri / 4 \ fool 

ChPT 0(p ) [3^ 
ChPT 0(/) [34] 
ChPT 0(/) [34] 
ChPT 0(/) [35] 


—4.9 ± 0.6 

-4.5 
-3.3 ±2.5 
-4.0 ±2.1 


5.2 ± 0.1 

5.9 
2.8 ± 1.1 
1.6 ±1.0 


2.1 
1.2 ± 1.7 


17 ± 10 

5.7 
3.5 ±0.6 


lAM 0{p'') [21] 


-3.7 ±0.2 


5.0 ±0.4 


0.8 ±3.8 


6.2 ±5.7 



TABLE IL Samples of other sets of LECs: First row: SU(3) 
analysis of nK scattering using Roy-Steiner equations. Sec- 
ond and third rows: K14 analysis to 0(p*) and 0{p^), respec- 
tively. Naively, we have combined quadratically the SU(3) 
LECs errors there. Fourth row: Roy Equations analysis. Un- 
certainties from imaginary parts and unknown O(p^) LECs 
combined quadratically. Last row, values used in [23] with the 
one-loop lAM. All LECs are evaluated at the scale /x = 770 
MeV 



these values as a crude indication of the size of systematic 
uncertainties. From the sets in [34j we note that, even 
for the same analysis, the values of the 0{p'^) LECs can 
be somewhat different whether they are obtained from 
a pure 0{p'^) calculation or including the 0{p^) correc- 
tions. Hence, it should not come as a surprise later that 
the 0{p^) values obtained from a unitarized fit, which 
includes part of the higher order corrections, may also 
come out somewhat different from the values obtained in 
a pure ChPT 0{p^) analysis. 

As a final comment concerning ChPT parameters, it 
is possible and usual to write the NNLO tttt scattering 
amplitude in terms of just six parameters 5i, . . . , fog; mul- 
tiplying each one of the energy dependent polynomials al- 
lowed by Lorentz invariance and chiral symmetry. Thus, 
the knowledge of 6 constants is enough to describe tttt 
scattering to that order. However, these hi parameters 
do carry a dependence on and the full knowledge of 
all the U and rj constants is needed to extrapolate to un- 
physical values of m,rj which is the object of this work, 
and the reason why we need to determine eleven param- 
eters instead of just six. 

Now, elastic unitarity implies for partial waves, at 
physical values of s, that: 



liat{s) = cr(s)lt(s)l^ 



rmlA(,s) = -a(s), (1) 



where ct{s) = 2p/-\/s, p being the center of mass momen- 
tum. As a consequence, the modulus of t{s) is related to 
its phase: 



t{s) ^ \t{s)\e'^^'^ ^ e'*'(^) sin,5(s)/a(s). 



(2) 



This "phase shift" (5(s), which determines completely the 
amplitude, is the usual way to parametrize partial waves, 
that we will use next to predict the amplitude variation 
when the pion mass is changed. Of course, before ex- 
trapolating to other pion masses, we will compare the 
ChPT amplitudes, with and without unitarization, with 
the existing experimental data. 

ChPT amplitudes, being an expansion, satisfy unitar- 
ity only pcrturbatively: 



Imt2=0, Imti ^ (j\t2\ , Imie = 2*2 R-et4 , 



(3) 
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In particular, ChPT partial waves are expected to violate 
unitarity as s increases, since they are basically polyno- 
mials in s. In section |V] we will use ChPT inside disper- 
sion relations to obtain amplitudes that, while respecting 
the ChPT expansion at low energies, satisfy unitarity and 
allow and provide a good description of experiment up to 
higher energies. 

After this brief introduction to ChPT and its notation, 
we are now ready to present our first calculations. 

III. RESULTS WITHIN STANDARD CHPT 

Using the equations above, the phase shift within stan- 
dard ChPT is obtained as a series expansion (see [36| for 
a prescription on how to perform this expansion): 

5 = cr{t2+RcU) + 0{p^), (4) 

2 

6 = c7{t2+RcU + Rcte) + -(ata)^ + O(p^), 

o 

which are the expressions used in our one-loop and two- 
loop calculations, respectively, that we detail next. 

Now, let us recall that the pion — and quark — mass de- 
pendence of the partial waves t{s) within ChPT comes 
from two different sources: from kinematics, through 
pion propagators, or from the dynamics encoded in the 
vertices. In particular, the threshold shift is purely of 
a kinematic nature and rather trivial to understand. 
Therefore, although tttt phase shifts are customarily pre- 
sented in terms of -y/s, we are showing them here as a 
function of the center of mass momentum p, which is 
also more convenient to compare to lattice studies. With 
this kinematic threshold effect "subtracted" , the remain- 
ing dependence is rather mild for most partial waves. 
As we will see this soft dependence of the Sij{p) on 
has been also found for 1 — 2 waves in very recent lattice 
calculations 



A. One-loop ChPT 

In Fig. [1] we show the phase shifts from the one-loop 
ChPT, i.e. 0(/), for the (/, J) = (0,0), (2,0), (1,1) tttt 
scattering waves. Note that for the (1,1) channel, the 
description fails much before p ~ 300 MeV. This mo- 
mentum is typically below the p{770) resonance region, 
which is a natural applicability bound for the ChPT se- 
ries. This resonance has a relatively narrow shape, cor- 
responding to a pole close to the real axis in the second 
Riemann sheet, which, of course, is completely missed by 
one-loop ChPT except in its very low energy tail. In con- 
trast, one-loop ChPT is giving a fairly good description 
of the (0, 0) channel even up to, say p — 350 or 400 MeV. 
In this case there is also a resonance — the scalar a (or 
/o(600)) — , but it is very wide and its pole is deep in the 
complex plane, so that it is not seen in the real axis as 
the typical sharp rise in the phase. For this reason, and 



despite being an expansion which has no such a pole in 
the complex plane, ChPT results are not very different 
qualitatively from the data in this channel. Finally, we 
see that the one-loop description of the (2, 0) channel is 
also reasonably good up to such high momentum, mostly 
due to the fact that this channel has no resonances and 
also that the data are not particularly precise. 

The gray areas in the figure cover the uncertainties due 
to the statistical error in the LECs detailed in the pre- 
vious section. In order to calculate these areas we have 
used a Montecarlo sampling. For each phase-shift calcu- 
lation we have generated 5000 different samples of LECs 
using a Gaussian distribution with variances equal to the 
errors quoted in Table HI To avoid a confusing overlap- 
ping between uncertainty bands, we only show the one 
corresponding to the physical pion mass. In the appendix 
we provide a detailed study of the evolution of these un- 
certainties with TOtt ■ As a general feature for both scalar 
and vector waves, the relative uncertainty of the phase 
at a given momentum grows slowly with m^r . 

Once we have checked where one-loop ChPT calcula- 
tions provide an acceptable description of data, we can 
now compare, also in Fig. [1] with the the results ob- 
tained if we change the pion mass from its physical value 
to m„ = 230,300 and 350 MeV. The first observation 
is that the sign of the phase derivative does not change 
when increasing the pion mass, at least up to 350 MeV, 
which means that the attractive or repulsive nature of 
each wave is conserved. 

In that figure we have represented with an arrow the 
direction of the phase movement as TOjt increases. Thus, 
the next observation is that both scalar phase shifts in- 
crease in absolute value as grows, whereas the phase 
of the vector channel decreases. 

The behavior of the phase at low momentum in the 
vector channel may seem surprising at first, because sev- 
eral lattice works [j-H, the chiral effective treatment j2^ , 
as well as the lAM [21|, predict that the p(770) mass in- 
creases much slower than the 27r threshold as grows. 
But then, when the yo(770) peak reaches a given momen- 
tum, the phase there should be 7r/2 to a very good ap- 
proximation. Therefore, one would expect naively the 
phase at low momentum to rise as grows. However 
the model independent ChPT analysis, tells us otherwise. 
We will see in detail in section ITVl why this intuitive pic- 
ture fails and the phase shift actually has to decrease at 
first and increase later on. 

Finally, in Fig. [5] we show the one-loop ChPT results 
for the D waves: (/, J) = (0,2) and (2,2). We show 
these separately because both them vanish at 0{p^), so 
that the one-loop 0{p'^) calculation is just their LO con- 
tribution. Actually they are both very small at low ener- 
gies. We can see in the figures that the one-loop ChPT 
calculation provides an acceptable solution for the (2,2) 
wave up to relatively high momentum, but obviously it 
cannot reproduce the resonance shape of the /2(1270) 
resonance in the (0, 2) channel. As before, we only show 
the uncertainty band due to the statistical errors on the 
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FIG. 1: S and P wave tttt phase shifts from standard ChPT 
up to one loop. Different hnes stand for different pion masses 
= 139.57, 230, 300 and 350 MeV, respectively. Since the 
lines are too close to each other, we only show error bands for 
the physical mass. Experimental data come from [37| (black 
circles) and the precise model independent dispersive data 
analysis from [3^ (white circles). The arrows show the direc- 
tion of increasing . 
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FIG. 2: D wave tttt phase shifts from standard ChPT up 
to one loop. Different lines stand for different pion masses 
= 139.57, 230, 300 and 350 MeV, respectively. Since the 
lines are too close to each other, we only show error bands for 
the physical mass. Experimental data come from [STj. The 
arrows show the direction of increasing mir. 



LECs for the physical pion mass, obtained again from a 
Montecarlo Gaussian sample. Relative uncertainties for 
different pion masses are detailed in the appendix. 

Note that, in contrast to the scalar waves, both tensor 
phase shifts decrease in absolute value as the pion mass 
increases not too far from its physical value. In this sense, 
they are more similar to the vector channel behavior. 
Remarkably, for larger pion masses and momentum the 
(0,2) phase shift even changes sign and the derivative 
becomes negative. However, this behavior is not found 
at two loops, as we will see in the next subsection. 



B. Two- loop ChPT 

We use the two-loop tttt scattering calculation in [27| . 
Note, however, that instead of the usual 6i . . . 6g param- 
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eters, in order to implement the ttitt dependence we need 
to use the one-loop li . . .I4 and parameters in Ta- 
ble m In Fig. [3] we show the resulting phase shifts for 
the (/, J) = (0, 0), (1,1) and (2,0) waves for the physical 
m^r but also for — 230, 300 and 350 MeV. 

The uncertainty bands, which we show only for the 
physical pion mass — see the appendix for other masses — 
are once again calculated with a Montecarlo Gaussian 
sampling of 5000 sets of LECs, using as standard devia- 
tions the uncertainties quoted in TablelH The only excep- 
tion are the ri...4 parameters, which are estimated from 
resonance saturation and, as in [29|, we have assumed 
that all values in the interval from to twice the estima- 
tion are equally likely. Of course, we want to emphasize 
that this is just an estimate of the values of the 0{p^) 
parameters, which are rather difficult to determine. Pos- 
sible improvements in their determinations could come 
from future lattice-QCD calculations, as it has already 
been done with the 0{p^) LECs (see [30| for a review) or 
from the use of recent dispesive data analysis like that in 
[ssj inside threshold parameter sum rules [3l|. 

Also, since the renormalization scale ^ where the esti- 
mates for ri...4 and Vf apply is not known, another source 
of uncertainty appears. Our calculations are made at 
ji — 770 MeV so, in order to account for the uncertainty 
due to that choice, we have followed [29] again and we 
have calculated the shift occurring in the phase shift if 
ri...4 are fixed and the scale is changed to ^ = 500 MeV 
and /i = 1 GeV. That shift is added in quadrature to the 
errors given by the Montecarlo sampling. 

The general features of the one-loop description still 
apply to the two-loop case. Namely, all waves keep their 
attractive or repulsive nature, and both scalar phases in- 
crease in absolute value as grows, whereas the vector 
channel phase decreases. The counterintuitive behavior 
of the p(770) is therefore a robust prediction of ChPT. 
In the next section we will explain with a simple model 
why chiral symmetry requires this behavior. Still, the 
description of the (0,0) wave is fair only up to p = 300 
or 350 MeV, although it has improved remarkably in the 
low energy region, where the data are most recent and 
reliable, as they come from Kg/^ decays. The (1, 1) phase 
is now much closer to the experimental data, and thus it 
seems to provide a fairly good representation up to, say 
p = 200 MeV. However, the description of the (2, 0) has 
deteriorated for higher momenta, and seems to be good 
only up to, roughly, 200 or 250 MeV. 

However, despite the qualitative m^r dependence being 
similar to the one-loop case, quantitatively the effect is 
stronger. In absolute value all phase shifts grow faster 
with to two loops than they did to one loop. 

In Fig. m we show the two- loop result for the D waves. 
As commented before, these waves have no 0{p^) term, 
so, this 0{p^) calculation is just a next to leading order 
calculation. We can see that the differences with the one 
loop case are dramatic. The (/, J) = (0, 2) phase suffers 
a remarkable improvement, being able to describe the tail 
of the 72(1275) resonance up to momentum of the order 
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FIG. 3: S and P wave tttt phase shifts from standard ChPT 
up to two loops. Different lines stand for different pion 
masses: continuous, long dashed, short dashed and dotted 
for Mn = 139.57, 230, 300 and 350 MeV, respectively. Since 
the lines are too close to each other, we only show error bands 
for the physical mass. Experimental data (black circles) come 
from [331 and the precise model independent dispersive data 
analysis from [s^] (white circles). The arrows show the direc- 
tion of increasing m,r . 
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FIG. 4: D wave tttt phase shifts from standard ChPT up 
to two loops. Different lines stand for different pion masses: 
continuous, long dashed, short dashed and dotted for m-^ — 
139.57, 230, 300 and 350 MeV, respectively. Since the lines 
are too close to each other, we only show error bands for 
the physical mass. Experimental data (black circles) come 
from 37]- The arrows show the direction of increasing m^. 



C. Comparison with lattice results for 1 — 2 and 

> 350 MeV 

As we have already commented, there are very recent 
lattice results on phase shifts for the / = 2, J = (l^.[l5j 
and J — 2 channels [14]. In Figs. [5] and IH] we compare 
the one and two-loop calculations within standard ChPT, 
first for the physical mass versus experimental data, and 
then for = 396, 420, 444 and 524 MeV, versus lattice 
results. 

When we examine Fig. [SJ corresponding to the 1 = 2, 
J = phase shifts, the first observation is that all 
lattice points with p < 200 MeV are well described 
within the uncertainties of one-loop ChPT, even up to 
= 444 MeV. From the figure, we observe that a 
pion mass of 524 MeV seems out of reach and will not 
be considered any longer. Beyond that momentum, the 
ChPT calculation bends downwards and misses all other 
lattice results with higher momenta. Remarkably the 
two- loop ChPT results do not improve this agreement. 
Actually, the two-loop calculation describes somewhat 
worse the lattice data and seems to move consistently 
to more negative values than those observed on the lat- 
tice, as rriTT grows higher. Let us remark that the cur- 
vature downwards is larger in the two-loop result than 
just to one loop. In view of the figures it seems that the 
standard ChPT applicability limit is, at best, somewhere 
around p ~ 150 — 200 MeV, up to of the order of 
400 - 440 MeV. 

Unfortunately, for the 1 = 2, J = 2 channel, shown in 
Fig.ini there are no lattice results available at low momen- 
tum. Surprisingly, the one-loop calculation agrees quite 
nicely with the lattice values up to around p ~ 500 MeV 
even for the highest pion mass. However, the two-loop 
results show a very strong to^ dependence that is in com- 
plete disagreement with the behavior predicted by the 
lattice simulations. Even the tendency is wrong, since 
the absolute value of the phase seems to grow with rriTr 
whereas lattice results may suggest a decrease. Let us 
nevertheless recall that for D-waves the tree level ampli- 
tude vanishes, so that one and two- loop calculations cor- 
respond only to leading and next to leading order results. 
Higher order calculations may be needed to improve and 
stabilize the D wave description. 



of 350 MeV. Contrary to the one-loop case, within the TOtt 
range of this study, the (0, 2) phase does not become neg- 
ative. Finally, the (2, 2) phase shift fails to describe even 
the sign of the data, and is only relatively close to the 
data points below 150 MeV. Furthermore, the one- loop 
rriTT phase shift dependence was opposite to the two-loop 
case: from more negative to less negative for the former 
versus from positive to negative for the second. The pre- 
dictions for this channel are therefore not very robust, 
which is also corroborated by the large uncertainties for 
higher m,^ that can be found in the appendix. 



IV. RECONCILING THE PHASE SHIFT AND 
RESONANCE BEHAVIOR IN THE VECTOR 
CHANNEL 

We have seen that, within ChPT, the low momentum 
phase shift of the vector channel is found to decrease 
as m,r grows. This is a model independent result and 
looks rather robust since it is obtained both at one and 
two-loops. However, lattice results 0-01, the chiral ef- 
fective treatment 



[25| . as well as the lAM [2lJ predict 
that, in terms of momentum, the p(770) peak gets closer 
and closer to threshold. Thus, for any low momentum 
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5: One and two loop standard ChPT phase shifts for 
= 2, J = channel compared to lattice results coming 
(circles) and [l5i] (triangles). 
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FIG. 6: One and two loop standard ChPT phase shifts for 
the 1 = 2, J = 2 channel compared to lattice results coming 
from [13] (circles). Note the large difference between one and 
two loop results. 



choice, and as ttItt increases, the p(770) peak reaches that 
given momentum, so that the phase there should be n/2. 
Therefore one would naively expect the phase shift for 
any fixed low momentum to grow with m^. 

Actually, this is what one would find if, to describe 
the p{770) resonance pole, one uses the very simple and 
intuitive (but, as we will see below, incomplete) Breit- 
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Wigner model 



t{s) 



-^sMT{p)/2p 
s - AP + iMT{p) 



where, = s/A — and the width is: 

3 



p_ 

PR 



(5) 



(6) 



where M is the resonance mass, pn is the pion momentum 
at the resonance energy pj^ — — so that t{s) 

behaves correctly at threshold, t{s) ^ p^K Note that 
Tr = g'^p^/enM'^ is the p(770) decay width. 

For the sake of simplicity, let us now assume that the 
resonance mass M and coupling remain constant when 
changing the pion mass uIt^. This implies that M and 
r(p) are independent. For our illustration purposes 
here, this is a fairly good approximation to what has been 
found on the lattice or with the lAM, and it could be 
considered as the leading order term in the expansion 
(see for the p(770) mass). 

In such case, however, the phase-shift m^r dependence 
near threshold does not follow what is obtained from 
ChPT (or the lAM as we will see below). In particu- 
lar, since 



MT{p) 
tan<5(p) = -7-2 



(7) 



the only m^r dependence in S (for a given p) is through 



PR (and d{pj^)/d{mi 



85 



dS 



-1) so that 
4Mr(p) 



d{p 



R.) 



(4p2 _ 4p|)^ -f NPT{py 



>0. (8) 



However, in ChPT, for low p we have shown in Figs. [T] 
and [3] that d5/d{ml) < 0. 

Of course, it is very well known that a simple Breit- 
Wigner vector formalism is not consistent with the chiral 
expansion unless there are some additional low energy 
contributions — or contact terms in the Lagrangian for- 
malism . Just to keep things very simple we can use a 
modification of the Breit- Wigner parametrization, which 
is widely used in analysis of tttt scattering and other phe- 
nomenology involving decays into light mesons [40| . and 
reads 



r(p) = r, 



PR/ 



Diipr) 



Di jPRr) 
Diipr) ■ 



(9) 



Here r(p) is m^r independent and Di(pr) are the Blatt- 
Weisskopf centrifugal barrier functions [4J] , that for I = 
1 read Di{pr) = 1 + {pr)^ . AU the TOtt dependence is 
carried by pR and the new parameter r, which is usually 
interpreted as a crude estimate of the "size" of the meson, 
although it should not be identified with its mean square 
charge radius. At low momentum we now find 



dSjp) ^ l+p%ir^y 
9(to^) 



^P% 



■Mfip), 



(10) 



where (r^)' stands for dr'^/d{m'^). In order to have a 
decreasing phase shift at low p when increasing m^, we 
just need 1 +PR{r'^y < 0. We will see below that this is 
actually required by chiral symmetry at leading order in 
the pion mass expansion. This would explain the phase 
decrease seen in ChPT for not too large ttItt, even though 
the p(770) is approaching threshold as m^r grows. Of 
course, when ttItt grows too large, and particularly in the 
limit when the p{770) tends to threshold, so that pr — 0, 
the derivative is positive, and the phase shift increases, 
as one would have expected naively. 

Let us then check that chiral symmetry actually re- 
quires 1 -|- p|j(r2)' < at least for low pion masses. We 
can estimate the leading dependence of r by compar- 
ing the low momentum and mass expansion of the am- 
plitude in Eqs.® using ©, with that of ChPT. In par- 
ticular, since in this simple model we have only one pa- 
rameter, r, we will only compare the scattering lengths. 
Our aim is just to reproduce the leading order TOtt depen- 
dence, since we have already made additional approxima- 
tions and simplifications (like the constancy of the /c(770) 
mass and coupling). We define the scattering length, a, 
as Re i ~ p^(a + bp^ + •••). The low p expansion of the 
amplitude in Eqs.(l5|) using ^ leads to 



m^MTR (1 + {pRrf 
4p%^ 



(11) 



Mp% 



(l + \M\' + 0{ml] 



This result has to be compared with that of ChPT: 
achPT = l/247r/2 + 0{ml). Matching with ChPT we 
obtain for P 



Pr 



1 



1 M 



O(m^) (12) 
O(m^) ~ (4.3 GeV"^)^ (13) 



The value obtained with this ChPT estimation is com- 
patible with what is found in the literature (r ^ 4 — 5 
GeV-i)) [40|. 

Note that the size r explodes as — > 0. However, 
this is a very well known feature of hadrons, at least 
for the charge radius. Actually, the squared charged ra- 
dius of the pion and the nucleon show a logm^ singu- 
larity m, Hll, Ell and the Pauli radius of the nucleon 
an additional l/m,r singularity [4l[. Nevertheless, as we 
have commented, our P parameter should not be directly 
identified with the p(770) charged radius, although our 
results suggest that they may have a similar singularity. 

With this ttItt dependence for r we find that 



l+P%ir'y^l 



Mp%, 



(14) 



which is negative for the physical values of the parame- 
ters. This guarantees that dS{p)/d{m'^) < for not 
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far from rnP^y , and sufficiently low p, as is obtained in 
ChPT. 

The decrease is a robust feature of ChPT, although the 
pure chiral expansion cannot reproduce the p(770) reso- 
nance. Of course, the model we have presented here is 
very simple and naive, but provides a qualitative and in- 
tuitive explanation of why chiral symmetry implies that 
the vector phase shift at low momenta first decreases, al- 
though it may increase later as ttItt grows. This model 
cannot be pushed too far. In particular, we cannot repro- 
duce the chiral behavior of the scattering length beyond 
leading order or even the slope parameter. 

It is however possible to incorporate simultaneously 
the p(770) pole and the full low energy ChPT expansion 
to one and two-loops. In the next section we will explain 
the technique in detail and later on we will show how 
it describes the existing lattice data up to much higher 
momentum than standard ChPT. Actually we will check 
how the vector phase shift decreases first and then in- 
creases as m^r grows. 



V. UNITARIZED CHPT: THE INVERSE 
AMPLITUDE METHOD 

As we have already commented in Sect. Ill Al the par- 
tial waves obtained from the ChPT expansion are ba- 
sically a truncated series in momenta or energies and 
cannot satisfy elastic unitarity, Eq.([T]), exactly, but only 
perturbatively, as in Eq.(|3]). 

There is, however, a well known technique, known as 
unitarization, to obtain expressions for partial waves that 
satisfy elastic unitarity, have the correct analytic struc- 
ture in terms of cuts in the complex plane, and simulta- 
neously respect the ChPT expansion up to a given order. 
Here we will make use of the elastic Inverse Amplitude 
Method (lAM) — or a slightly modified version — that im- 
plements the fully renormalized one or two-loop ChPT 
expansion at low energies but does not introduce any 
spurious parameter in the unitarization procedure. Had 
we used other, possibly simpler but very successful, uni- 
tarization techniques with spurious parameters, like cut- 
off or any other regulator, we should have had to worry 
about the unknown dependence of that scale. 

The lAM [l^ uses elastic unitarity and the ChPT ex- 
pansion to evaluate a once subtracted dispersion relation 
for the inverse amplitude. The analytic structure of 1/t 
consists on a right cut from threshold to oo, a left cut 
from — oo to 0, and possible poles coming from zeros of t. 
We can write then a once subtracted dispersion relation 
for 1/t, the subtraction point being sa, 

t{s) IT Jjic {s' - sa){s' ~ s) (15) 
+ LC{l/t) + PC{l/t), 

where LC{l/t) stands for a similar integral over the left 
cut and PC{l/t) is the contribution of the pole at sa- 



The choice of sa is, in principle, arbitrary, but since we 
want to use the information encoded in the ChPT se- 
ries, we are then limited to the low energy region, prefer- 
ably, below threshold. Now, scalar waves vanish at the so 
called Adler zero that lies in the real axis below thresh- 
old and in practice this is a very convenient choice for 
Sa, which has actually motivated our notation. For other 
waves there is no such an Adler zero, and the subtraction 
point can be taken, for instance at s = 0. It is important 
to remark that the choice of subtraction point, as long 
as it lies between the left and right cut, has only a very 
small numerical effect [l^ on the physical region. Up 
to here everything is exact. The most relevant observa- 
tion is that, following Eq. ([1]), on the elastic cut we know 
exactly Iml/t = —a. 

Now we are going to derive the lAM within one-loop 
ChPT. First, the Adler zero position can be approxi- 
mated as, Syi = S2 + S4 + • • • , where ^2 vanishes at S2, 
ti -V i4 vanishes at s^ + S4, and so on. On the right cut 
we can evaluate exactly Iml/i = —a = — Imi4/f2, as 
can be read from Eqs. ([T]) and Since the left cut is 
weighted at low energies we can use one-loop ChPT to 
approximate LC{l/t) ~ LC(— ^4/^2)- The pole contribu- 
tion PC{l/t) can be safely calculated with ChPT since 
it involves derivatives of t evaluated at sa, which is a 
low energy point where ChPT is perfectly justified. Al- 
together ^we arrive to a modified one-loop lAM (mlAM) 
formula [ig |: 

,mIAM _ ^2 

.mlAM , , X (g2 - Sa){s - S2)[4(S2) - ^4(^2)] 

A ^ti{S2) , 

S-SA 

(16) 

where the prime denotes the first derivative with respect 
to s and where we use for sa in the numerical calcula- 
tions its NLO approximation S2 -I- S4. The standard lAM 
formula is recovered for A"^^^^'^ = 0, which is indeed the 
case for all partial waves excegt the scalar ones. In the 
original lAM derivation [l3, [3 A"^^^^'^ was neglected 
since it formally yields a higher order contribution and is 
numerically very small except near the Adler zero. How- 
ever, if ^"J-f-^M neglected, the lAM Adler zero occurs 
at S2, correctly only to LO, is a double zero instead of a 
simple one, and a spurious pole appears close to the Adler 
zero. All of these caveats disappear with the ml AM, and 
the differences between the lAM and the mlAM in the 
physical and resonance region are of the order of 1%. 

It is important to remark that ChPT has not been used 
at all for calculations of t{s) for positive energies above 
threshold. Note that the use of ChPT is well justified to 
calculate sa, and PC{l/t), since these are low energy cal- 
culations. ChPT has also been used to calculate the left 
cut integral, which, despite extending to infinity, is heav- 
ily weighted at low energies, which once again justifies 
the use of ChPT. The left cut and the elastic approxi- 
mation are the only approximations used to obtain the 
lAM, but no other model dependent assumptions have 
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been made. In particular there are no spurious parame- 
ters included in the lAM derivation, but just the ChPT 

LECs, TOjr and f-rr- 

Remarkably, these simple equations (either the lAM 
or the mlAM) ensure elastic unitarity, match ChPT at 
low energies and, using LECs compatible with existing 
determinations, describe fairly well data up to somewhat 
less than 1 GeV, generating the p, K* , a and k resonances 
as poles on the second Riemann sheet [l8|. 

The extension to two loops is very similar and straight- 
forward for the lAM or the mlAM [M]: 



%IAM 



to 



+ te{s2) + 
-4(52) + 



2U{s2)t'i{s2) <1(S2) 



^2(^2) t'^{s2){s - S2) 

(S - S2){SA - S2) 



SA 



t'2{s2)~Us2) 



t\{s2? +t'i{s2)U{s2) 
t'2is2) 



(17) 

Let us now remark that both in the one and two-loop 
derivations above, we have assumed that ^2 is not iden- 
tically zero. However, this is only the case for scalar and 
vector partial waves. Unfortunately, as seen in Eq.(j31), 
when ^2(3) = the first imaginary part appears at 0{p^), 
namely, at three loops. Therefore we cannot recast the 
dispersion relation in terms of the full ChPT expansion 
unless we make use of ts{s), a calculation that does not 
exist. In (43 ), and using only the term of the form cs*, 
it was shown that the /2(1275) shape could be fairly well 
fitted with the lAM and a c value of the correct order 
of magnitude expected from dimensional grounds. This 
was justified because the /2 resonance appears at high 
s >> and the other 0{p^) terms, containing pion 
mass powers, could be neglected. However, in this work 
we want to make m,r much larger than its physical value 
and we need the m^r dependence. It is therefore not so 
well justified to neglect all the terms except cs^. For 
that reason we are limited to use the lAM for scalar and 
vector partial waves. 

Hence, using the lAM or the mlAM, we can study how 
the generated p and a poles evolve by changing to^t in 
the one loop lAM amplitudes [2l| or two-loop amplitudes 
[23 |. and describe the dependence of their masses, widths 
and couplings on m^r. In [2l| the mlAM was used for 
the p and a chiral extrapolation, because, for the scalar 
and at high m^, one resonance pole gets near the lAM 
spurious pole, a problem that is nicely solved with the 
mlAM. Nevertheless, in the physical region and near the 
other generated poles, the differences between lAM and 
mlAM approaches are almost negligible, even for high 
pion masses. 

Of course, the poles are not the only object of study on 
the lattice. Actually, lattice results are already available 
for phase shifts in / = 2 channels, where no pole exists. 



Moreover these channels were not studied in [2l|, |2J] . It 
is also very likely that lattice results on phase shifts for 
other channels will be available soon. For these reasons 
we will now let m^r vary within our unitarized ChPT 
expressions, with the aim of extending the phase shift 
predictions based on ChPT, up to higher masses and mo- 
menta. 



VI. RESULTS WITH THE lAM AND CHPT 

Let us first recall, as already explained in some of the 
very first works on the I AM [18l|, and repeated in many 
other instances [l^, IH, El], that when the central val- 
ues of the standard LECs are used, the lAM only im- 
proves ChPT up to a couple of hundred MeV higher 
and resonances are only reproduced qualitatively. For 
a semi-quantitative description of resonances, which is 
what we will do next, one has to fit the data and the 
resulting LECs are slightly modified from those obtained 
from pure ChPT. Since the lAM contains contributions 
that count as higher order in ChPT (in particular the nu- 
merically relevant s-channel logarithms) , one would very 
naively expect the LECs from the one-loop lAM to lie 
somewhere in between the one and two loop values from 
ChPT. This is actually observed, since the 0{p*) lAM 
LECs in Table IIIII lie somewhere between the one and 
two loop analysis of pure ChPT listed in Table [TTl al- 
though closer to the ChPT O(p^) analysis in the two 
first rows of that table. In contrast, the 0{p'^) values of 
the LECs for the two-loop lAM in Table |IV] are closer 
to the two loop analyses like that in Table |T] or those 
in the third and fourth row of Table |TT1 Let us empha- 
size that the variation between the 0{p'^) LECs values 
between the one and two loop analyses already occurs in 
pure ChPT — particularly for Zj. The I AM simply follows 
a similar pattern. 

Before changing the pion mass, let us note that for 
the lAM we are assuming the elastic approximation and 
therefore, when increasing ttt-tt we should allow for some 
TTTT elastic regime, which is guaranteed if ttXtt < 500 MeV, 
although it has been found that relatively stable unita- 
rized results can be obtained for all waves only up to 
TOjr ~ 300 — 350 MeV [13] . Of course some waves are 
more stable than others. In particular the elastic lAM 
approximation is quite good up to larger energies for the 
(i, J) = (2, 0) (roughly up to ^/s ~ 1200 - 1300, see 
p5|). since it has no resonances and does not couple to 
KK. We will actually check that for this channel we can 
stretch the applicability range and still get fairly good 
agreement with recent lattice results for relatively large 
pion masses. 



A. One-loop lAM 

In Fig.[7]we show the lAM results to one-loop in ChPT, 
using the LECs in Table IIIII obtained by an updated fit 
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0[P ) 






-3.9 ±0.2 




4.3 ±0.4 


I3M 


0.18 ± 1.11 




6.17 ± 1.39 



TABLE III: LECs used in this work for the one-loop lAM, 
obtained from a fit to the dispersive data analysis of ^3q. Both 
^3 and I4, axe fixed to the standard values given in Tabled The 
scale is set to ^ = 770 MeV. 



to the output from the recent and precise dispersive data 
analysis in [38], and fixing ^3 and Z4 to the updated val- 
ues in Table m The uncertainties are mostly systematic, 
arising from different choices of the maximum energy up 
to where we make the fit of the (0, 0) channel, which 
we have chosen between 500 and 800 MeV, the other 
channels are fitted up to 1 GeV. Note that the resulting 
LECs are consistent within one standard deviation with 
the results we used in [2l|, that we list in the last row 
of Table HT] We first note that the experimental data 
is fairly well described up to the region where inelastic 
effects (or resonances like the /o(980)) become relevant. 
This includes the /9(770) resonance shape, but also the 
wide shape of the /o(600). The gray bands in the fig- 
ures cover the uncertainties in our results obtained from 
a Monte Carlo Gaussian sampling of the U statistical er- 
ror bars also listed in the table. As usual, and to avoid 
confusion due to many overlapping gray bands, we only 
show the uncertainty for the physical pion mass. Details 
on uncertainties for higher masses can be found in the 
appendix. 

The general features for the scalar isoscalar channel 
are very similar to the one-loop non-unitarized results. 
Namely, the phase shift conserve its positive sign and 
increases in absolute value as TOjt grows. 

However, the 1 = 2 channel behavior is rather differ- 
ent. First, the dependence is even milder than for 
the non-unitarized case. In the very low momentum re- 
gion, roughly below p = 200 MeV, the phase increases 
in absolute value as it happened with standard one-loop 
ChPT. However, for larger momentum, the m^r depen- 
dence is the opposite, and the phase starts decreasing its 
absolute value. As we will see later on, this is the be- 
havior found on recent lattice results, which cannot be 
reproduced by a crude extrapolation of one-loop ChPT 
to larger momentum. 

Something similar occurs in the vector channel, al- 
though enhanced by the presence of the p(770) resonance 
that ChPT failed to reproduce. Now we see that the 
phase increases as the two pion threshold grows and gets 
closer to the resonance. This is the intuitive behavior one 
would expect when getting close to the resonance. How- 
ever, one should observe that it is not incompatible with 
the phase decrease observed in standard ChPT at low 
energies. To see this, in Fig. [HI we show a blow up of the 
very low energy region of the vector channel, where we 




p (MeV) 

FIG. 7: S and P wave tttt phase shifts from unitarized 
ChPT up to one loop. Different lines stand for different pion 
masses: continuous, long dashed, short dashed and dotted for 
Mn = 139.57, 230, 300 and 350 MeV, respectively. Since the 
lines are too close to each other, we only show error bands for 
the physical mass. Experimental data come from ^7} (black 
circles) and the precise model independent dispersive data 
analysis from [s^l (white circles). The arrows show the direc- 
tion of increasing mir. See Fig. [8] for a blow up of the low 
momentum region of the 7 = 1, J = 1 phase shift. 
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FIG. 8: TTTT / = 1, J = 1 phase shift from unitarized 
ChPT up to one loop. The continuous line stands for = 
139.57 MeV and the dotted line for = 350 MeV. Sim- 
ilarly to the ChPT case, in the low momentum region the 
phase shift decreases as grows. However at higher mo- 
mentum it increases with the pion mass, due to the presence 
of the p(770) resonance. 



can see that the JAM behaves similarly to ChPT, namely, 
the phase decreases as ttItt grows. As explained before, 
this only happens in the very low momentum regime, 
since, as seen in the figure, for higher momentum the 
phase shift increases again since the JAM is able to re- 
construct the p(770) resonance, which is closer and closer 
to threshold as grows. 

In the next subsection we will see that these gen- 
eral features and improvements with respect to non- 
unitarized ChPT are even more dramatic when consider- 
ing the two-loop calculation. 



B. Two- loop I AM 

In Figs. 1^ and [TUl we show the results of the two-loop 
JAM for the two best fits in [11], "A" and "D" , whose cor- 
responding sets of LECs we provide in Table HV] These 
fits have been obtained from an lAM fit to experimental 
data but also to lattice results on f-^ , Mp and the isospin 
2 scattering length. Note that by fitting only the experi- 
mental data one determines better the LECs that govern 
the s dependence, but not so well those governing the 
771^ dependence. That is the reason why some existing 
lattice results on Mp and the 1 — 2 scalar scattering 
length were also included in the fits of 24]. Unfortu- 
nately, the experimental data in the resonance region are 
frequently in conflict with one another, and to a lesser 
extent, something similar happens for the lattice results 
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TABLE IV: Low energy constants obtained from a fits [2J| to 
experimental data on elastic tttt scattering and lattice results 
on /tt, Mp and the isospin 2 scattering length as well as a l/7Vc 
leading behavior of a pure qq state. Many of these sets are not 
quite compatible with each other and suffer large systematic 
uncertainties. These two fits correspond to different ways 
of weighting the existing experimental and lattice data sets, 
which are detailed in [2J]. The values correspond to the scale 
II = 770 MeV. 



mentioned above. Fits A and D correspond to different 
ways of weighting the conflicting experimental and lattice 
results, including some educated estimates for systematic 
uncertainties. The details can be found in [23|. These flts 
give rather stable results for all observables in the elas- 
tic region, up to — 300 — 350 MeV, and somewhat 
beyond for some particular waves, like (/, J) = (2, 0). 

Note that the qualitative behavior of all waves is sim- 
ilar in Figs. [9] and [TOl The difference between flt A and 
D is purely quantitative: in flt A the dependence is 
just stronger than in fit D. 

Remarkably, almost all the features described for the 
one-loop unitarized case remain in the two-loop unita- 
rized flts. Quantitatively there are small differences, 
since the dependence at two loops seems somewhat 
stronger in the scalar waves, and somewhat weaker in 
the vector channel. This somewhat stronger m^r depen- 
dence produces the only signiflcant, and relevant, differ- 
ence with the one-loop lAM. Both the one and two-loop 
lAM generate the /o (600) or a resonance as a pole deep in 
the complex plane, which mass grows much slower than 
the two pion threshold, so that the '"bump"' that this 
wide resonance produces in the (0, 0) phase is bigger and 
gets closer to threshold. Actually, as shown in j21|| the 
two conjugated poles of the /o(600) move in the second, 
unphysical, Riemann sheet, until they reach the real axis 
below threshold, where the two poles are no longer conju- 
gated. As TO^ keeps growing one of them jumps into the 
first Riemann sheet below threshold becoming a bound 
state. By Levinson's theorem [i^l, this implies that the 
phase at threshold increases by tt. For the lAM to one- 
loop this jump occurs for larger than 350 MeV, but 
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FIG. 9: S and P wave tttt phase shifts from the two-loop lAM 
"fit A" in . The conventions are as in Fig. [7] The arrows 
show the direction of increasing m,r. The difference between 
these curves and those in Fig llUl are an indication of the order 
of magnitude of our uncertainties. 



FIG. 10: S and P wave tttt phase shifts from the two-loop 
lAM "fit D" in [11]. The conventions are as in Fig. [T] The 
arrows show the direction of increasing m-n- The difference 
between these curves and those in Fig|5]are an indication of 
the order of magnitude of our uncertainties. 
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since the 171^^ dependence is stronger for the lAM at two 
loops, this jump can aheady be seen in Figs . 1^ and [TUl for 
the ttItt = 350 MeV, which behavior thus reflects the ex- 
istence of a bound state. Let us emphasize that the same 
behavior would be observed to one loop — although for 
higher — but it will never be seen in standard ChPT, 
which cannot generate a pole. 

However, when comparing with the non-unitarized 
two-loop results in Fig. [31 we see that unitarization not 
only improves the vector channel by describing the p(770) 
resonance, but also the 1 — 2 channel is nicely described 
up to much higher momentum, even though this channel 
is non-resonant. We will profit from this lack of com- 
plicated resonant structures in the I ~ 2 scalar wave, 
and also from the fact that this channel does not couple 
to KK^ to extrapolate to higher pion masses where we 
will see that the unitarized results are in much better 
agreement than standard ChPT with some recent lattice 
results. 



VII. COMPARISON WITH LATTICE RESULTS 

FOR / = 2 AND m„ > 350 MEV 

In FigdU we show the results from the one and two- 
loop JAM with very recent results on the lattice [13, [ll| 
for the 1 = 2 scalar channel. Note that the data below 
p = 200 MeV is still fairly well described by the lAM, as 
it happened with ChPT, but that the lAM is not bending 
down and getting away from higher momentum data as 
it happened with standard ChPT results. Actually, the 
lAM results follow qualitatively the shape of the lattice 
data. Moreover, the ttItt dependence is much milder than 
for plain ChPT, in better agreement with the findings on 
the lattice. Let us remark that we do not aim at precision 
here because pion masses of 400 MeV are probably close 
to the JAM applicability bound. Our approach should 
become more reliable below 300-350 MeV, where we ex- 
pect lattice results to appear soon. Still, the remarkable 
improvement with respect to the standard ChPT results 
is pretty clear. 

As previously commented, the JAM cannot be directly 
applied to the D waves, since their tree level contribution 
vanishes. Further modifications of the lAM would be 
needed which are beyond the scope of this work. 



VIII. SUMMARY AND DISCUSSION 

In this work we have studied the pion mass dependence 
of TTTT elastic scattering phase shifts. 

On the one hand we have presented results for one and 
two-loop standard Chiral Perturbation Theory using a 
set of LECs obtained from a dispersive analysis in the 
literature. We have seen that this first approach is, of 
course, limited to low momentum, say below 300 MeV, 
depending on the channel, and pion masses up to 400- 
450 MeV. For the scalar and vector waves, we have found 
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a rather stable behavior between the one and two loop 
calculations within that momentum range. We have seen 
that at this very low momentum, the absolute value of 
scalar phase shifts increases as the pion mass grows, so 
that these channels enhance their attractive or repulsive 
nature. We have found that up to momenta less than 
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200 MeV, the ChPT results are in fair agreement with 
lattice data for the scalar 1 = 2 channel. 

We have found that, surprisingly, the vector phase- 
shift at very low momentum decreases as TOtt grows 
within the applicability region. This may seem coun- 
terintuitive, since from lattice and other effective theory 
techniques, as m^r grows one expects the two-pion thresh- 
old to approach fast the p(770) mass. We have neverthe- 
less shown with a very simple and intuitive model why 
very basic requirements about chiral symmetry impose 
such a decrease on the phase for low momentum and not 
too large TOtt. 

We have also shown results within standard ChPT for 
the angular momentum 2 phase shifts. These are much 
less stable when comparing one and two loop results. 
Particularly for the (/, J) = (2, 2) channel, the one and 
two loop results show an opposite behavior, and the two 
loop calculation is also at odds with the TOtt dependence 
found on the lattice. Of course, one has to keep in mind 
that for D waves the one and two-loop calculations corre- 
spond to leading and next to leading order calculations, 
contrary to scalar and vector channels, where they cor- 
respond to next to leading and next to next to leading 
calculations. It is very likely that higher order calcula- 
tions, or better determinations of LECs, which are highly 
correlated, may improve this situation for D waves. 

Finally, we have used ChPT inside a dispersion rela- 
tion to extend the analysis of scalar and vector waves to 
higher momentum by means of the so called Inverse Am- 
plitude Method. This unitarization technique describes 
remarkably well the data up to energies of the order of 
1 or 1.2 GeV, depending on the channel and has been 
shown to describe well the TOtt dependence of several ob- 
servables like Mp, or the 1 = 2 scalar scattering length. 

The description provided by this method is of course 
compatible with that of standard ChPT at very low mo- 
mentum. However at higher momentum it reconstructs 
the behavior of the p{770) resonance, which, for a given 
choice of low momentum, translates into a decreasing 
phase for smaller but a growing phase for larger 
until the p(770) mass coincides with that particular mo- 
mentum choice. In addition, we have shown that the 
unitarized 1 = 2 scalar phase shift has the correct quali- 
tative behavior for momentum beyond 200 MeV. Despite 
being close to the applicability bounds of the approach, 
we have actually shown that the lAM beyond p = 150- 
200 MeV improves dramatically the description of lattice 
results with respect to ChPT and explains their very mild 
TO^ dependence. 

Intuitively, the phase shift evolution of the SO and P 
channels is dominated by the presence of the /o(600) and 
p(770) resonances and their pion mass dependence, stud- 
ied in |21, 23 25]. Since the masses of both resonances 
seem to grow slower than the pion mass, they come closer 
and closer to threshold, so that, naively one would expect 
the interaction to grow stronger and the phase to raise 
once the resonance is sufficiently close to the momentum 
where the phase is measured. Actually, this is what is 



found for the SO channel, whose phase raises noticeably 
as ttItt grows. At the limit of the range of applicability 
of the two- loop lAM, the /o(600) even becomes a bound 
state and by Levinson's theorem we see the phase to in- 
crease by TT at threshold. However, the naive expecta- 
tions may not be met if the resonance is still not close 
enough to threshold. In such case, the phase may seem 
to decrease at first due to the finite size of the resonance, 
which effect has been illustrated in a simple model of 
the p{770). Only when the p{77Q) is sufficiently close 
to threshold, the naively expected behavior is observed. 
Concerning the S2 wave, we have found a very mild m,r 
dependence for the phase shift, when expressed in terms 
of the momentum, in good agreement with recent lattice 
calculations. This can be understood from the absence 
of resonant structures in this channel. Of course, ChPT 
can only reproduce the low energy tails of the resonances, 
which we have generated by means of ChPT unitarized 
with the I AM. For the D waves, the I AM cannot be ap- 
plied to this order, and we have to rely on ChPT only. 
However, the behavior observed can also be understood 
from the presence of the /2(1270) resonance in the DO 
channel, and a similar behavior to the p(770) in its own 
channel. For the D2 channel the ChPT results are not 
sufficiently precise to make any conclusive statement. 

Apart from understanding the dependence of these ob- 
servables on QCD parameters on the pion mass, we con- 
sider that this work is of interest as a guideline for future 
studies of lattice QCD. 
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Appendix A: Phase shift uncertainties for different 

In Fig. [T^ we plot the relative uncertainties of the stan- 
dard ChPT phase shift calculation. As we have already 
seen, standard ChPT is limited to low momentum and 
thus we only show momentum up to p = 300 MeV. For 
the scalar and vector waves we see that in the low mo- 
mentum region the errors grow with the pion mass. This 
is in agreement with the fact that the LECs that gov- 
ern the mass dependence of the partial waves carry the 
biggest uncertainties. For D-waves, the relative uncer- 
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tainty is much bigger than for lower angular momentum 
waves. (Note the difference in scales between the D waves 
and the rest of the plots). This is due to the fact that for 
D-wavcs the tree level calculation vanishes and therefore 
the one and two loop results are just leading and next 
to leading order. In the case of S02 to one loop the error 
seems to explode for the highest masses due to the phase 
shift changing from a positive to a negative value in the 
region of interest. The same occurs for S22 to one loop 



for the physical value of the pion mass. Finally, the value 
of 622 to two loops changes from negative to positive for 
the lightest masses of the pion. 

In Fig. [13] we show the relative uncertainties for the 
lAM phase shifts. We find again that for scalar waves 
they grow bigger as the pion mass is increased. The same 
happens for the vector phase shift below the jo(770) peak. 
The highest uncertainty on Sn occurs when the slope of 
the phase shift reaches its maximum value. 
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FIG. 13: TTTv phase shift errors normahzed to the value of the phase shifts in unitarized ChPT to one loop. Different hnes 
stand for different pion masses: continuous, long dashed, short dashed and dotted for = 139.57, 230, 300 and 350 MeV, 
respectively. 
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